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Strawberry cultivar identification based on hypervariable SSR markers 
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We genotyped strawberry cultivars by two newly selected and two previously reported SSR markers. All four 
markers produced interpretable electropherograms from 75 accessions consisting of 72 Fragaria x ananassa 
cultivars or lines and three octoploid Fragaria species accessions. These SSR markers were highly polymor- 
phic; in particular, one of the newly developed markers, FxaHGA02P13, was capable of distinguishing all 
of the accessions except for a mutant strain that was derived from another accession in the set. When two 
markers were combined, all 48 full-sib individuals could be distinguished. Fingerprinting patterns were 
reproducible between multiple samples, including the leaves, sepals, and fruit flesh of the same accession. 
Principal-coordinate analysis of the 75 accessions detected several groups, which reflect taxon and breeding 
site. Together with other available markers, these SSR markers will contribute to the management of straw- 
berry genetic resources and the protection of breeders' rights. 
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Introduction 

Cultivated strawberry, Fragaria x ananassa (In = 8x = 56), 
is an economically important crop worldwide. In Japan, the 
value of strawberry production was more than 168 billion 
yen in 2007, making strawberry the third most valuable crop 
after rice and tomato (Ministry of Agriculture, Forestry and 
Fisheries 2010; http://www.maff.go jp/j/tokei/index.html). 

Strawberry cultivars are usually vegetatively propagated 
from runners. Recently, the unregulated propagation and 
distribution of patented cultivars has become a serious prob- 
lem (Kunihisa 2010). To protect breeders' rights and to 
maintain the genetic integrity of genetic resources, several 
molecular marker systems have been developed for straw- 
berry cultivar identification. Kunihisa et al. (2003, 2005, 
2009) developed 25 cleavage amplified polymorphic se- 
quence (CAPS) markers and confirmed their ability to dis- 
tinguish among 117 cultivars, except for mutant strains 
derived from others in the set. This system is highly 
reproducible, produces results that are simple to interpret, 
and does not require an autosequencer for the detection of 
polymorphisms. Meanwhile, because the degree of poly- 
morphism of each marker was low, a number of markers 
were required for cultivar identification. In addition, CAPS 
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marker analysis requires the use of restriction enzyme diges- 
tions. For these reasons, analyses of large numbers of sam- 
ples can be time-consuming. Tasaki et al. (2008) developed 
three sets of multiplex PCR primers based on RAPD (ran- 
dom amplified polymorphic DNA)-STS (sequence tagged 
site) markers. A combination of the three sets could dis- 
tinguish 25 cultivars, though the number of cultivars dis- 
tinguished by any individual primer set was only 1 to 4. 
Because the authors analyzed only 25 cultivars, it was 
unknown whether the marker system could distinguish addi- 
tional cultivars. 

Recently, the number of strawberry cultivars registered in 
Japan has substantially increased. Whereas only 35 cultivars 
were registered from 1980 to 1989, 53 were registered from 
1990 to 1999, 120 were registered from 2000 to March 2011, 
and applications for registration have been published for 20 
more. To ensure the sensitive and rapid identification of the 
increasing number of cultivars, we expect that additional 
polymorphic markers will be also useful. 

In general, simple sequence repeats (SSRs) are highly 
polymorphic, reproducible, and easy to assay; for these rea- 
sons, they have been widely used for cultivar identification 
in clonally propagated crops such as grape, apple, pear, po- 
tato, carnation, turfgrass species, and others (Ashkenazi et 
al. 2001, Bowers et al. 1996, Guilford et al. 1997, Kimura et 
al. 2002, 2009, Smulders et al. 2003, Wang et al. 2010, 
Yamamoto et al. 2006). Also for strawberry, development of 
SSR markers suitable for cultivar identification is a current 



Strawberry cultivar identification by SSR markers 



421 



area of interest. Shimomura and Hirashima (2006) devel- 
oped four SSR markers for strawberry, but the degree of 
polymorphism of those markers was relatively low; even 
when the most polymorphic marker was used, 2 of the 10 
cultivars tested showed the same genotype and could not be 
distinguished. Govan et al. (2008) selected 10 SSR markers 
suitable for the fingerprinting of strawberry cultivars, and 
Brunings et al. (2010) tested the applicability of these 10 
markers for the genotyping of Florida strawberry cultivars. 
Brunings et al. (2010) found that two markers, EMFvl04 
and EMFvil36, were both capable of distinguishing all ac- 
cessions when combined with one additional marker. In 
these studies (Brunings et al. 2010, Govan et al. 2008), 
'Nyoho' was the only Japanese cultivar analyzed, so the ap- 
plicability of those markers for discriminating Japanese cul- 
tivars was not determined. 

In this study, we revealed the fingerprinting patterns of 
75 Fragaria accessions, including a number of Japanese cul- 
tivars, by using two newly developed SSR markers (T. 
Nunome et al, manuscript in preparation) and two SSR 
markers (EMFvl04 and EMFvil36) reported in previous 
studies (Brunings et al. 2010, Govan et al. 2008). We also 
evaluated the discrimination power of these markers for 
highly related individuals by analyzing a set of 48 full-sibs. 
By using the detected SSR polymorphisms among the 75 
accessions, we analyzed the genetic relationships within this 
group. This study is one of the first to identify highly 
polymorphic SSR markers suitable for the identification of 
Japanese strawberry cultivars. 

Materials and Methods 

We collected leaves from 75 accessions consisting of 72 
F. x ananassa cultivars or lines and three octoploid Fragaria 
species accessions (Table 1), which were grown at the 
National Agricultural Research Center for Tohoku Region. 
For most of the accessions, we collected leaves from multi- 
ple plants to confirm reproducibility. For four cultivars 
('Miyazaki-natsuharuka', 'Natsuakari', 'Summer-berry' and 
'Summer-princess'), we also collected sepal and fruit flesh 
samples to confirm the reproducibility of the genotype when 
different tissue samples were used. To estimate the discrim- 
inating power of these SSR markers for highly related indi- 
viduals, we also collected leaves from 48 full-sib individuals 
produced from a cross between 'Sagahonoka' and 'Summer- 
berry'. Total DNA was extracted by using a modified PEG 
method (Rowland and Nguyen 1993) with Plant DNAzol 
Reagent (Invitrogen, Carlsbad, CA, USA), as described by 
Sugimoto et al. (2005), or by using a DNeasy Plant Mini Kit 
(Qiagen, Hilden, Germany). 

Each of the samples was genotyped by using four SSR 
primer pairs (Table 2). The primers for EMFvl04 and 
EMFvil36 were previously reported by Govan et al. (2008). 
FxaHGA02P13 and FxaAGA21Fl 1 were selected from 200 
primer pairs, which were newly developed by Nunome et al. 
(in preparation), according to the ease of PCR amplification 



and interpretation, reproducibility and degree of polymor- 
phism. For each sample, the polymerase chain reaction 
(PCR) was performed by two methods to examine the ro- 
bustness of the fingerprint patterns: the original method of 
Schuelke (2000) and a slightly modified method based on 
the method of Schuelke (2000) (hereafter, "modified meth- 
od"). In the original method (Schuelke 2000), the PCR reac- 
tion mix contained 5.0 uE of GoTaq Colorless Master Mix 
(Promega, Madison, WI, USA), 1.6pmol of the reverse 
primer, 1.6pmol of the universal fluorescent (D2, D3, or 
D4)-labeled M13(-21) primer, 0.4 pmol of the forward prim- 
er with a M13(-21) tail at its 5' end, and 10 ng template DNA 
in a total volume of 10 uE. A leader sequence (GTTTCTT) 
was appended to the 5' end of the reverse primers to mini- 
mize the appearance of stutter bands in the electrophero- 
grams (Brownstein et al. 1996). Thermocycling conditions 
were as follows: 5 min at 94°C; then 30 cycles of 30 s at 
94°C, 45 s at 56°C and 45 s at 72°C; followed by 8 cycles of 
30 s at 94°C, 45 s at 53°C and 45 s at 72°C and a final ex- 
tension at 72°C for 10 min. In the modified method, the 
fluorescent-labeled primer was added after the first PCR am- 
plification, which was performed in a 10-uL volume. The 
amplification was performed as follows: 5 min at 94°C, fol- 
lowed by 10 cycles of 45 s at 94°C, 1 min at 52°C and 45 s at 
72°C. Then, 0.8 uL of 2pmol/uE of the fluorescent-labeled 
primer was added to each PCR tube. The second amplifi- 
cation was performed as follows: 15 s at 94°C; 30 cycles 
of 45 s at 94°C, 1 min at 45°C and 45 s at 72°C and a final 
extension step at 72°C for 7 min. These PCR reactions 
were carried out in a GeneAmp PCR System Model 9700 
(Applied Biosystems, Foster City, CA, USA) or a PCR 
Thermal Cycler Dice (Takara, Tokyo, Japan). The PCR 
products were run on a CEQ 8000 autosequencer (Beckman- 
Coulter, Fullerton, CA, USA). 

We analyzed the relationships between accessions based 
on the genotypes of the four SSRs in the following way. For 
each accession, any peak detected in at least one accession 
was scored as present (1) or absent (0) across all of the SSR 
markers. Then, the binary data were combined into a single 
matrix and used for the calculation of a genetic similarity 
matrix with the aid of PAUP 4.0 (Swofford 2002). Based on 
the matrix, we conducted a principal-coordinate (PCO) anal- 
ysis to visualize the relationships among the accessions by 
using the software PCO 1.0 (Iwata 2004). For marker 
EMFvil36, we used the genotype data obtained by the mod- 
ified method, because more peaks could be detected (see 
Results and Discussion). 

Results and Discussion 

In this study, we genotyped strawberry cultivars by two new- 
ly selected SSR markers and two previously reported SSRs 
(Brunings et al. 2010, Govan et al. 2008). All four markers 
revealed interpretable electropherograms from all acces- 
sions surveyed. The genotype of each accession is listed in 
Supplemental Table 1 . When we performed genotyping by 
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Table 1. Name and country of origin of strawberry accessions analyzed in this study 



Accession No. 


Accession name 


Country 


Accession No. 


Accession name 


Country 


1 


Aiko 


USA 


39 


Morioka 32 


Japan 


2 


Akihime 


Japan 


40 


Morioka 33 


Japan 


3 


Akitaberry 


Japan 


41 


Morioka 34 


Japan 


4 


Amaotome 


Japan 


42 


Morioka 35 


Japan 


5 


Benihoppe 


Japan 


43 


Natsuakari 


Japan 


6 


Beruruju 


Japan 


44 


Nohime 


Japan 


7 


Blakemore 


USA 


45 


Nyoho 


Japan 


8 


Bolero 


UK 


46 


Ohishi-shikinari 


Japan 


9 


Cambridge Favourite 


UK 


47 


Otomegokoro 


Japan 


10 


Cardinal 


USA 


48 


Oze-akarin 


Japan 


11 


Dekoruju 


Japan 


49 


Pajaro 


USA 


12 


Deutch Evern 


Germany 


50 


Pechka 


Japan 


13 


Donner 


USA 


51 


Pegasus 


UK 


14 


Elsanta 


Netherlands 


52 


Pelican 


USA 


15 


Everberry 


Japan 


53 


Raiho 


Japan 


16 


Florence 


UK 


54 


Red Gauntlet 


UK 


17 


Fukuba 


Japan 


55 


Reiko 


Japan 


18 


Harumi 


Japan 


56 


Sachinoka 


Japan 


19 


Haruyoi 


Japan 


57 


Sagahonoka 


Japan 


20 


Hatsukuni 


Japan 


58 


Sequoia 


USA 


21 


Hecker 


USA 


59 


Summer-berry 


Japan 


22 


Hitachihime 


Japan 


60 


Summer-candy 


Japan 


23 


Hogyoku 


Japan 


61 


Summer-drop 


Japan 


24 


Hokowase 


Japan 


62 


Summer-fairy 


Japan 


25 


HS 138 


Japan 


63 


Summer-princess 


Japan 


26 


Ibarakiss 


Japan 


64 


Summer-tiara 


Japan 


27 


Kaorino 


Japan 


65 


Tioga 


USA 


28 


Karenberry 


Japan 


66 


Tochihime 


Japan 


29 


Kitaekubo 


Japan 


67 


Tochihitomi 


Japan 


30 


Kitanokagayaki 


Japan 


68 


Tochiotome 


Japan 


31 


Komachiberry 


Japan 


69 


Toyonoka 


Japan 


32 


Mae 


UK 


70 


Tsuburoman 


Japan 


33 


Marshall 


USA 


71 


Uzushio 


Japan 


34 


Miyazaki-natsuharuka 


Japan 


72 


Yayoihime 


Japan 


35 


Miyoshi 


Japan 


73 


Fragaria chiloensis 'PI551445' 




36 


Moikko 


Japan 


74 


F. virginiana #1" 




37 


Morioka 16 


Japan 


75 


F. virginiana #2 b 




38 


Morioka 30 


Japan 









" Seeds were purchased from B&T World Seeds, Paguignan, France, in 2002. 

* Germplasm conserved in the National Agricultural Research Center for Tohoku Region, Morioka, Japan. 



the original PCR amplification method of Schuelke (2000), 
the number of peaks scored per marker ranged from 20 
(FxaAGA21Fll) to 34 (EMFvl04) and the number of 
unique genotypes ranged from 61 (FxaAGA21Fl 1) to 74 
(FxaHGA02P13; Table 2). When we performed genotyping 
by the modified method, two additional peaks (145 and 155) 
were detected at EMFvil36. The fingerprinting patterns de- 
tected by FxaHGA02P13, FxaAGA21Fll and EMFvl04, 
did not change, suggesting that these three SSR markers are 
more robust and less susceptible to the difference in PCR 
method than EMFvil36. 

One of the newly developed markers, FxaHGA02P13, 
was capable of distinguishing 74 unique genotypes among 
75 accessions (Table 2). Thus, FxaHGA02P13 has high dis- 
criminating ability comparable to that of EMFvl04, which 



showed the highest discriminating ability in previous studies 
(Brunings et al. 2010, Govan et al. 2008). 'Akitaberry', a so- 
matic mutant that originated from 'Morioka 16', had the 
identical genotypes to that of 'Morioka 16' across all four 
SSR markers, as was seen in the previous CAPS analysis 
(Kunihisa et al. 2009). Although the number of unique geno- 
types detected by FxaAGA21Fll (61 among 75 acces- 
sions, i.e., 81.3%) was lowest among the four markers sur- 
veyed here, the discriminating rates for 10 markers selected 
from over 100 SSRs by Govan et al. (2008) ranged from 
83.3% to 100.0% in their study, and were much lower in 
Florida cultivars (Brunings et al. 2010), indicating that 
FxaAGA21Fll has discriminating power comparable to 
those previously selected markers. Although a comprehen- 
sive analysis using multiple markers is more authoritative, 
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Table 2. The four SSR primer pairs used in the analysis of 75 strawberry accessions, the number of peaks scored, the range of observed peak 
sizes, and the number of unique genotypes detected 



SSR name 


Forward primer sequence 


Reverse primer sequence 


No. of peaks 
scored 


Range of 
observed peak 
sizes (bp) 


No. of unique 
genotypes 
detected 


FxaHGA02P13 


CCAGGCGCTTGGTCTTGTACTACT 


CCCATTTCCCCCAAATCTAACAAT 


25 


244-302 


74 


FxaAGA21Fll 


CAATTCACAATGGCTGATGACGAT 


GCACTCAGACATATTTTGGGAGGG 


20 


134-185 


61 


EMFvl04 


TGGAAACATTCTTACATAGCCAAA 


CAGACGAGTCCTTCATGTGC 


34 


98-166 


73 


EFMvil36 


GAGCCTGCTACGCTTTTCTATG 


CCTCTGATTCGATGATTTGCT 


29 (27)° 


140-197 


65 



0 Two peaks at EMFvil36 were not detected when PCR amplification was performed by the original method of Schuelke (2000). All other geno- 
types were the same with both methods. 



analysis with even a single SSR marker can provide a fair 
amount of information. 

Although peak sizes increased corresponding to the 
lengths of the sequences we added to the primers, the geno- 
types at EMFvl04 and EMFvil36 obtained by the modified 
method almost corresponded to those obtained in the previ- 
ous study (Govan et al. 2008) for cultivars analyzed in com- 
mon. However, some discrepancies such as appearance or 
disappearance of peaks and slight peak size changes were 
observed, as were also seen in Brunings et al. (2010). Tech- 
nical variations between studies were considered to be one 
of the major reasons of such discrepancies. Meanwhile, 
highly reproducible patterns were observed among experi- 
mental methods used in a previous study (Govan et al. 
2008). Also in this study, reproducible amplification pat- 
terns were observed, including between multiple tissue types 
of the same accession (Fig. 1), and between two PCR meth- 
ods for three markers as mentioned above. Therefore, al- 
though it is required to pay careful attention to the reproduc- 
ibility of genotype, especially in the comparison of 
genotypes obtained by different methods or laboratories, 
useful and convenient genotyping can be performed by 
choice of appropriate experimental method and SSR mark- 
ers. For data harmonization between methods or laboratories 
(Doveri et al. 2008, Govan et al. 2008, This et al. 2004), 
genotypes of several standard cultivars, and several fre- 
quently observed peaks (e.g., 244 in FxaHGA02P13, 146 
and 150 in FxaAGA21Fll, 102 in EMFvl04 and 149 in 
EMFvil36) will serve as landmarks. 

Further, we tested the ability of these four markers to dis- 
criminate among 48 full-sib individuals produced from a 
cross between 'Sagahonoka' and 'Summer-berry', to esti- 
mate the resolving power of these markers for highly related 
individuals. In the set of 48 full-sibs, the number of unique 
genotypes detected by each marker ranged from 4 
(EMFvil36) to 35 (EMFvl04; Table 3). These values had 
the same relative ranking as the number of peaks that segre- 
gated among the full-sib progeny: which ranged from 3 
(EMFvil36) to 9 (EMFvl04). When two markers were com- 
bined, the discrimination power increased greatly: by com- 
bining any two of the three markers FxaHGA02P13, 
FxaAGA21Fl 1 and EMFvl04, all 48 full-sib individuals 
could be distinguished. Thus, although it is necessary to pay 



'Summer-berry' 




'Summer-princess' 




240 250 260 270 280 290 300 310 (bp) 



Fig. 1. Amplified fragment patterns of SSR FxaHGA02P13 from dif- 
ferent tissue samples of two cultivars. The number above each peak in- 
dicates the fragment size. 

attention to the degree of differentiation among parents, the 
combination of just two markers enabled us to discriminate 
between even highly related individuals. 

In addition to cultivar identification, hypervariable SSR 
markers may be also useful for the elucidation of the genome 
composition of the octoploid strawberry. Two main hypoth- 
eses for the genome composition of F. x ananassa have been 
proposed: (1) AAAA'BBB'B' (Bringhurst 1990, Kunihisa 
et al. 2005) and (2) AAA'A'BBBB (Senanayake and 
Bringhurst, 1967, Lerceteau-Kohler et al. 2003). Highly 
polymorphic SSR markers, which are theoretically possible 
to amplify eight products from four homoeologous loci, may 
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Table 3. The number of unique genotypes detected in the analysis of 
48 full-sibs by single or pairs of SSR markers used in this study 



FxaHGA02P13 


FxaAGA21Fl 1 


EMFvl04 


EMFvil36 


FxaHGA02P13 18° 


48 


48 


43 


FxaAGA21Fl 1 


28" 


48 


45 


EMFvl04 




35" 


47 


EMFvil36 






4° 



" Value for the single SSR marker. 

afford direct evidence for mode of inheritance in the cultivat- 
ed strawberry. For example, in the genotyping of offspring 
of a cross (AjA 2 A 3 A 4 A 5 A 6 A 7 A 8 x Z\Z\ Z{L\ ZtZi Z{L\), 
either one of the two peaks of a parent (e.g., Ai and A? and 
so on) certainly inherited in the offspring but did not appear 
simultaneously in an offspring, these two peaks were in- 
ferred to be alleles in a pair. If such four pairs (A1A2/A3A4/ 
A 5 A 6 /A7A 8 ) are obtained, it could be one of the direct evi- 
dence for the hypothesis AAA A'BBB'B'. In this study, two 
peaks of EMFvl04 were considered to be alleles in a pair 
from the analysis of above-mentioned 48 full-sib individu- 
als, but unfortunately, other pairs were not identified due to 
the overlap of the same size peaks between parents. By se- 
lecting mating parents whose peaks are not overlapped, alle- 
les in a pair may be more clearly identified. 

PCO analysis of the 75 accessions detected several 
groups, which reflect taxon, country of origin, and breeding 
site in Japan, in which limited breeding materials have been 
repeatedly used to achieve flowering habit suited to the re- 
gion (Fig. 2). Japanese early-flowering June-bearers were 
plotted toward the positive end of the PCOl axis, which ac- 
counted for 11.6% of the total variation. Along the PC02 
axis, which accounted for 7.3% of the variation, groups corre- 
sponding to wild Fragaria species accessions, cultivars bred 
in America or Europe, and Japanese cultivars could be iden- 
tified. In Japan, strawberry breeding programs have been 
conducted to produce cultivars adapted to each region or 
cropping type. In the warmer regions of Japan, mainly early- 
flowering June-bearers have been selected because they are 
adapted to forcing culture, in which fruits are usually har- 
vested beginning in the early winter. On the other hand, in 
the cooler regions of Japan, mainly late-flowering June- 
bearers have been selected because they are better adapted to 
semi-forcing or open culture, in which fruits are harvested in 
spring and early-summer. In cooler regions, the breeding and 
use of everbearing cultivars is also popular. As a result of the 
repeated use of breeding materials suited to particular breed- 
ing objectives, differences in allelic composition between 
the early-flowering June-bearers and the late-flowering 
June-bearers and everbearers may have occurred. Several 
everbearers and late-flowering June-bearers that were de- 
rived from early- flowering June-bearers, such as 'Summer- 
princess' (accession 63; offspring of early- flowering 
'Nyoho'), 'Raiho' (accession 53; offspring of a seedling of 
'Nyoho') and 'Oze-akarin' (accession 48; a grandchild of 
early-flowering 'Tochiotome'), were plotted in the proximi- 



q O.OR+Ofl 




-3.0E-01 

-3.0E-01 0.0E+00 3.0E-01 

PCOl 

□ American and European cultivars 

O Japanese cultivars: early-flowering June-bearers 

A Japanese cultivars: late-flowering Junc-bcarcrs 

X Japanese cultivars: everbearers 

X Wild octoploid Fragaria species 

Fig. 2. Principal-coordinate (PCO) analysis of 75 accessions consist- 
ing of 72 F. x ananassa cultivars or lines and three octoploid Fragaria 
species accessions analyzed by four SSR markers. The first (PCOl) 
and second (PC02) axes explain 11.6% and 7.3% of the total varia- 
tion, respectively. Numbers correspond to accession numbers in 
Table 1. 



ty of early-flowering June-bearers. 

Three accessions of wild octoploid species F. chiloensis 
and F. virginiana (accessions 73, 74 and 75) were plotted at 
the edge of the group in Fig. 2. Although F. x ananassa is an 
interspecific hybrid of F. chiloensis and F. virginiana, only 
a few genotypes of these species contributed to the establish- 
ment of F. x ananassa (Dale and Sjulin 1990). Probably for 
this reason, the F. chiloensis and F. virginiana accessions 
used in this study showed a different allelic composition 
from the F. x ananassa accessions. Old cultivars such as the 
first Japanese cultivar, 'Fukuba' (accession 17; released in 
1899), 'Marshall' (accession 33, released in 1893), which 
was the oldest among cultivars analyzed in this study, 
'Cambridge Favourite' (accession 9; released in 1947) and 
'Blakemore' (accession 7; released in 1929), were plotted 
closer to wild Fragaria accessions than almost all of the 
other more recent cultivars. Although further investigation 
with more markers is needed, such pattern may reflect the 
historical transition of genetic composition in strawberry 
cultivars. Also, the fact that most of the early-flowering June- 
bearing cultivars in Japan after 'Toyonoka'(accession 69; 
registered in 1984) and 'Nyoho' (accession 45; registered in 
1985) converged in the right side of Fig. 2, while older early- 
flowering June-bearers such as 'Fukuba', 'Hogyoku' (acces- 
sion 23), 'Hatsukuni' (accession 20), 'Uzushio' (accession 
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71), 'Haruyoi' (accession 19), 'Hokowase' (accession 24) 
and 'Reiko' (accession 55) were scattered around the posi- 
tion, may suggest the genetic homogenization in recent 
years. 

This study identified highly polymorphic SSR markers 
which can be used for the identification of Japanese straw- 
berry cultivars. Together with other available markers, these 
SSRs will contribute to the management of genetic resources 
of strawberries and the protection of breeders' rights. 
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